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Abstract 

In this paper we determine, in dimension three, the effective conductivities of non peri- 
p ^ odic high-contrast two-phase cylindrical composites, placed in a constant magnetic field, with- 

out any assumption on the geometry of their cross sections. Our method, in the spirit of the 
H-convergence of Murat-Tartar, is based on a compactness result and the cylindrical nature of the 
i^h , microstructure. The homogenized laws we obtain extend those of the periodic fibre-reinforcing 

case of [17] to the case of periodic and non periodic composites with more general transversal 
geometries. 

1 Introduction 

> 
CN 

{S) ' At the end of the 19th century, it was discovered [2 1] that a constant magnetic field h modifies the 
symmetric conductivity matrix a of a conductor into a non symmetric matrix a(h). This is know as 
the Hall effect. In the Maclaurin series of the perturbed resistivity (<r(/i)) -1 the zeroth-order term 
coincides with the resistivity cr~ l in the absence of a magnetic field [ ]. In dimension two, h is 
a scalar and the first-order term is an antisymmetric matrix proportional to hJ; the coefficient of 
proportionality is called the Hall coefficient. In dimension three, /i£l 3 and the first-order term, in 
the Maclaurin series of (cr(/i)) -1 , is of the form S'(Rh) where <2>(£)j '■= £ x j and R is a 3 x 3 matrix 
called the Hall matrix [16]. In this work, we consider the idealized situation when the induced non 
symmetric part is proportional to the applied magnetic field: cr(/i) = a/3 + /3<c?(/i), where a and (3 
are two constant real numbers. For a given sequence of perturbed conductivities a n (h), it is of 
great interest, in electrodynamics [27, 32], to understand the influence of the magnetic field h on the 
effective Hall coefficient or the effective Hall matrix through the homogenization of cr n (h). 

Let us first review a few of the mathematical theory of homogenization of elliptic partial differ- 
ential equations of the form 

f -div(<T n Vn n ) = / in 
1 u n = on <9Q, 

where £1 is a bounded open subset of M 3 , a n is a sequence of matrix-valued functions in L°°(fi) 3x3 
and / is an element of i2" _1 (fl). This topic has been intensively studied for the last four decades 
providing a wide literature [34, 30, 31, 2]. In the context of conduction, when the conductiv- 
ity matrices a n are uniformly bounded, Spagnolo [34] with the G-convergence theory, Murat and 
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Tartar [30, 31] with the .ff-convergence theory showed that the solution u n 6 Hq(Q) of the conduc- 
tivity problem (1.1) strongly converges in L 2 (Q), up to a subsequence of n, to the solution of a limit 
conductivity problem of the same nature. The case of high-contrast conductivities is very different 
since non classical phenomena, such as nonlocal terms, may appear in the limit problem as shown, 
for instance, in [19, 25, 1, 18, 11, 26]. This does not happen in dimension two if the sequence a n 
is uniformly bounded from below. Briane [ ] and Casado-Dfas & Briane [13] proved that in that 
case the class of equations (1.1) is always compact in the sense that the limit equation of (1.1) is 
always of the same type. In [13] they proved some extensions of the well-known div-curl lemma of 
Murat-Tartar [ ] and deduce several compactness results under the assumption of equicoerciveness 
coupled with the L -boundedness of the sequence of conductivities. 

In this paper we are interested in the homogenization of a class of three-dimensional conductivity 
problems of the type 

| -div (a n (h)X7u n ) = f in O, 
I u n = on dQ, 

where a n (h) is an equi-coercive sequence of high-contrast two-phase conductivities perturbed by a 
constant magnetic field h £ M 3 of the form a n (h) := (1 — lQ ri )o"i(/i) + ^n„^2,n(h) where a2, n (h) is 
the unbounded perturbed conductivity of the highly conducting phase Q n and <7\{K) is the perturbed 
conductivity of the phase surrounding tt n . 

In dimension two, for the case of low magnetic field, Bergman [3] was the first author who came 
up with a general formula for the effective Hall coefficient of a periodic composite material in terms of 
the local Hall coefficients and some local currents solving the conductivity equations in the absence 
of a magnetic field. We refer also to the works [28, 4, 15] for other two-dimensional composites, 
to [5, 7, 8, 21, 22] for composites with microstructure independent of one coordinate (the so-called 
columnar composites) and to [ , 9] for the case of strong magnetic field. 

Recently, in dimension two, M. Briane and the second author [17] obtained the effective perturbed 
conductivity a*{h) of a sequence of isotropic high-contrast two-phase conductivities a n (h) in the case 
of strong magnetic field, i.e., when the symmetric part and the antisymmetric part of the conductivity 
are of the same order. By extending a duality principle from [14] and using a suitable Dykhne 
transformation, which (following Milton [28, 29]) changes non symmetric matrices into symmetric 
ones, they proved that the symmetric part of the effective perturbed conductivity er*(/t) is given in 
terms of the effective conductivity in the absence of a magnetic field. They subsequently compared 
their two-dimensional results to a three-dimensional periodic one and showed that the way a magnetic 
field perturbs the conductivity of a composite depends on the dimension. In order to compute the 
explicit perturbation formula in dimension three, they restricted themselves to a particular periodic 
fibre-reinforced structure, i.e., a structure completely described by any two-dimensional cross section 
transversal to the fibres (first introduced by Fenchenko, Khruslov [19] to derive a non local effect 
in homogenization). To our knowledge, only few results are known on the homogenization of both 
high-contrast and non symmetric conductivities in dimension three. 

The aim of this paper is to determine the effective perturbed conductivity of (1.2) for non 
periodic high-contrast two-phase cylindrical composites without any assumption on the geometry of 
the transversal microstructure. 

We first investigate the periodic case, that is, when a n (h)(-) = Y> n {h){- / e n ) where £ n (/i)(-) is 
a Y-periodic matrix-valued function and e n — > represents the size of the heterogeneities in the 
composite. In order to avoid non local effects in the limit problem, following Briane [12], we assume 
the existence of a sequence of positive numbers c n such that e\ c n tends to zero, as n goes to infinity, 
and satisfying the weighted Poincare-Wirtinger inequality 

2 

VVeH l (Y), J \E n (h)(y)\ V- f V dy dy < c n ! E n {y)VV ■ VVdy. 
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For a fixed n£N*, using the theory of exact relations of Grabovsky, Milton, Sage [21, 20] (thanks to 
the independence of the microstructure of the variable X3), we obtain the H-limit (<7 n ) associated 
with the periodic homogenization [2] of the oscillating sequence H n (-/e) as e — > 0. Then, we show 
that the sequence of constant conductivities ((7 n ) t converges to some a*(h) which, according to [12], 
coincides with the homogenized conductivity associated with the limit problem of (1.2). The obtained 
effective conductivity a*{h) is explicitly computed in terms of the homogenized conductivity <7*(/i) 
of the conduction problem posed in the (x±, X2)-plane transversal to the columnar composite (see 
Proposition 2.1). 

Most of the arguments and tools used in the periodic case crucially lie on the periodic nature 
of the microstructure. Therefore, a fundamentally different approach is necessary for the analysis 
of (1.2) when a n (h) is not periodic. 

In order to study the asymptotic behavior of the problem (1.2) in the non periodic case, using a 
method, in the spirit of the H-convergence of Murat- Tartar, we determine the limit, in an appropriate 
sense, of the current a n (h)\7u n . The key ingredient of this approach is a fundamental compactness 
result (see Lemma 3.1) based on a control of high conductivities in thin structures through weighted 
Poincare-Wirtinger type inequalities. This compactness lemma, combined with the two-dimensional 
results of [17] and the cylindrical structure of the composite allows us to obtain an explicit formula 
of a*(h), once again, in terms of the transversal homogenized conductivity a*{h) and of some bounded 
function 6 which, in some sense, takes account of the distribution of the highly conducting phase Vl n 
in Vt (see Theorem 3.1). 

The structure of the paper is the following: In Section 1.1 we set up some general notations. 
Section 2 deals with the periodic case. In Section 3 we extend the periodic result of Section 2 to a 
non periodic framework. Section 4 is devoted to some examples illustrating both the periodic and 
non periodic perturbation formulas. 

Here, we give some general notations and definitions. 

1.1 General notations and definitions 

• Q is a bounded open subset of M 3 with a Lipschitz boundary. The unit cube (—5, ^) 3 of M 3 
is denoted by Y. 

• For any subset co of $7, we denote by uJ the closure of u in R 3 . 

• e n is a sequence of positive real numbers converging to zero as n goes infinity. 

• For any matrix a in ]R rfx<i , a T denotes the transpose of the matrix a while a s denotes its 
symmetric part. For any invertible matrix a in M rfx<i , o"~ T := (c^ 1 ) = (c T ) 

• Id denotes the unit matrix in ~R dxd and J := (i q)- 

• For any h £ M 3 , <§(K) denotes the 3x3 antisymmetric matrix defined by $(h)x := h X x, 
for x € M 3 . 

• For any a,rj £ M rfxd , a < r] means that for any £ £ R rf , <r£ • £ < rjl; ■ £. 

• For any vector £ £ M 3 , £ £ M. 2 denotes the vector of its first two components 

e:=(6,6) T 

• V- denotes the gradient operator in M 3 with respect to the three variables (x\, X2, £3) while V- 
is the gradient operator in M 2 with respect to the first two variables (xi, £2): for any u £ H 1 ^), 
the function Vtt is defined on Q by 



where Q is the projection of f2 on the (x%, X2)-plane. 
• For any 3x3 matrix a, we denote by a the 2x2 matrix defined by 



Til T12 

a :■■ 

CT21 T22 



• The scalar product of two vectors u and v of R d is denoted by u ■ v. 

• | • | denotes, the euclidean norm in M. d , the subordinate norm in M. dxd and the Lebesgue 
measure. 

• For a Borel subset and a function u E L 1 ^) the average value of u over ui is denoted 
by 

it dx := - — - / u dx. 

When uj = Y, we simply denote this average value by (•). 

• We denote by 1^ the characteristic function of the set u. 

• We denote by ^ c (f2) the set of continuous functions with compact support in Vl. The subspace 
of ^c(O) of infinitely differentiable functions with compact support in VL is denoted by 

• We denote by ^q(O) the space of continuous functions on $7 vanishing on the boundary c?f2 
of Q endowed with the usual norm. 

• For any locally compact subset X of R , A4(X) denotes the set of Radon measures defined 
on X. 

• A sequence in is said to weakly-* converge to a measure [i if 

/ ipfi n (dx) ► / (p/j,(dx), for any ip G ^(ft). 

Jn n ^°° Jn 

• The space of Y-periodic functions which belong to L? (M ) (resp. H{ oc (R d )) is denoted 
by L%{Y) (resp. H^(Y)). 

• o(8) denotes a term of the form 8(^(8) where the limit of ((8) is zero, as 8 goes to zero. For any 
sequences (on)ngN* &nd (&n)n6N*; a n ~ means that a n = b n 

• Throughout the paper, the letter c denotes a positive constant the value of which is not given 
explicitly and may vary from line to line. 

In the sequel, we will use the following extension of H-convergence for two-dimensional high- 
contrast conductivities introduced in [13] for the symmetric case and extended in [14] to the non 
symmetric case: 

Definition 1.1. Let O be a bounded domain of M 2 and let a n £ L°°(S7) 2x2 be a sequence of equi- 
coercive matrix-valued functions. The sequence a n is said to H(A4(Q) 2 ) -converge to a matrix-valued 
function a* if for any distribution g in i/~ 1 (f2) 7 the solution u n of the problem 

{div (<7„,Vu n ) = g in £7, 
u n = on dQ, 



satisfies the convergences 



in Hlip), 



0"nVn n — cr*Vn weakly-* in Ai(Jiy 
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where u is the solution of the problem 



J div (a* Vit) = g in Q, 
1 u = on dfl. 

Let Q be a bounded open subset of M. 2 with a Lipschitz boundary and U n be a sequence of open 
subsets of f2. Let f2 be the bounded open cylinder £7 := x (0, 1) and f2 n the sequence of open 
cylinders fi n := f2 n x (0, 1). Consider a\ > 0, /?i G 1R and two sequences a 2) n > ot-x and /32, n G K. 
Define, for any li£l 3 , the two-phase isotropic conductivity 

( a\{h) := aih + P\£(h) in \ U n , / o -h 3 h 2 

o~ n (h) := < where S(h) : = I h 3 o -ft 



<?2,n{ h ) := a 2,nh + fo,n£{h) in ft n , \ ->»2 &l 







In the domain £1, the matrix- valued function o~ n (h) does not depend on the variable X3 and model the 
conductivity of a columnar heterogeneous medium. The phase fi n is the one of high conductivity: 02,™ 
and /32, n ar e unbounded. In order to ensure the L 1 (il) 3x3 -boundedness of the conductivity, we assume 
that the volume fraction of the highly conducting phase 9 n := |Q| _1 |il n | converges to zero and that 
the convergences 



dn(X2,n > «2 > 0, 

n— )-oo 

0n#2,n > ft£R, 



;i.3) 



hold. Assumption (1.3) can be rewritten 

n (T2,n( h ) = na 2 ,nh + O n fo, n £(h) > (7 2 (h) := Q 2 / 3 + fa£(h). 

Our aim is to study the homogenization of the Dirichlet problem, for / G i? _1 (17), 

f -div (<x n (h)Vun) = f in 0, 

\ u n = on an. { ' ' 

On the one hand, we consider the case of a periodic cylindrical composite without any assumption 
on the geometry of its cross section. This framework extends the one of the three-dimensional result 
of [17] where the highly conducting zone is a set of circular fibres. On the other hand, by the 
means of a compactness result (see Lemma 3.1), we analyse the case of cylindrical but non periodic 
composites. In both cases, we impose conditions, adapting [12], that prevent from the appearance 
of non local terms so that the limit equation of (1.4) is a conductivity one. 

In the sequel, we will omit the dependence on h of G\(h), o~2 t n(h) an d (72(h) denoting sim- 
ply CJi, cr 2 ,n and cr 2 . 



2 The periodic case 

In this section, we study the influence of a constant magnetic field h G M 3 on the effective conduc- 
tivity of a composite material where the highly conducting phase is periodically distributed but, con- 
trary to [ 7], the cross section of which has a general geometry. Consider a sequence u n = uj n x (0, 1) 
where u n is a sequence of subsets of (0, l) 2 with \u> n \ converging to 0, as n tends to infinity. Let Q n 
be the sequence of open subsets of f2 defined by 

0„ = tt n (J e n (u n + k) . 
The conductivity of the heterogeneous medium occupying Q is given by 

a n (h)(x) = Z n (h) , Viefi, (2.1) 
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where S n (/i)(-) is a ^-periodic function defined by 

E n (h) = a f J i + b n *(h) with { ?: =t Z 1 ^ Un +?' n } Un ' (2-2) 

For a fixed n G N*, let {tJ n )*{h) be the constant matrix defined by 

VAeR 3 , (<t„)*(/i)A = (£„(/i)W n A ), (2.3) 
where, for any A G M 3 , W^ 1 is the unique solution in H^(Y) of the auxiliary problem 

div (Z n (h)VW£) = in @'(R 3 ) and (W* - X ■ y) = 0, (2.4) 
which is equivalent to the variational cell problem 

f (E„(fc)vw* • v$) = o, v*eflf(y), 

1 (^ A (y)-A-y} = 0. 



(2.5) 



The matrix (c n )*(/t) is the homogenized conductivity of the oscillating sequence S n (-/e) as e — > (see, 
for instance, [2] for more details). 

The limit problem of the high-contrast three-dimensional equation (1.4) where o~ n (h) is given 
by (2.1) may include non local effects. In order to avoid such effects, we assume, following [12], that 
the weighted Poincare-Wirtinger inequality 

WGi^OO, Ja n V-jV < C nj a n \VV\ 2 , (2.6) 

holds true with 

e 2 n C n > 0. (2.7) 

n— >oo 

Under the assumptions (2.6) and (2.7), it was shown in [12] that the sequence of problems (1.4) 
converges to a conduction one with a homogenized conductivity cx*(/i). 

The main contribution of Proposition 2.1 below is to provide a formula for the effective conduc- 
tivity a*{h) of a cylindrical periodic composite the cross section of which has a general geometry. 

Proposition 2.1. Consider the sequence of problems (1.4) where a n (h) is the conductivity defined 
by (2.1)-(2.2). Assume that (1.3), (2.6) and (2.7) are satisfied. Then, there exists a constant ma- 
trix a*(h) such that, up to a subsequence, the solution u n of (1.4) weakly converges in Hq(CI) to the 
solution u of 

- div (er*(/i)Vw) = / in CI, 
u = on dCl. 



(2.8) 



Moreover, the homogenized matrix a*(h) is the limit of (er n ) (h) (see (2.3) ^ and is given by 

a*(h):=( a * T P *\, (2.9) 




where 



p* 

q* 

a* 



-[I3ih + h{<y*-ai)a 2 l ] Jh, 

[Pih + foa^^-a^fjh, (2.10) 
a\ + a 2 + /?2 ? 2 _1 (5 : i + <?2 — cr*)?^ 1 Jh ■ Jh, 



and, for any i = 1,2, 




Remark 2.1. For the sake of simplicity, throughout the paper, the symmetric part of c n (h) is sup- 
posed to be isotropic. However, the results we obtain can be extended to composites the components 
of which have anisotropic conductivities. 

Remark 2.2. It was shown in [12] that, due to the L 1 (y) 3x3 -boundedness of the se- 

quence is bounded. Thanks to (2.6) and (2.7), Theorem 2.1 of [12] ensures that the 

limit <7*(/i) obtained in the following way 



V n (h)(-) — > ><r,(h), 

\£J £->0 n-i-oo 

satisfies the convergence 



a n (h)Vu n — >■ a*(h)Vu in .M(f2) 3 , 
and then, coincides with the homogenized conductivity matrix in the problem (2.8). 

Remark 2.3. Since Q n has a columnar structure, the sequence <r n (h) given by 

<xi(/t 3 ) = ail 2 + PxhzJ infi\fi n , 



a n (h) := a n (h 3 ) 



depends only on the transversal variable (xi,X2) and is then associated with the two-dimensional 
problems, for any g £ 

J -div (a n (h 3 )Vv n ) = g in O, 
[ v n = on dil. 

Similarly to (2.3), we define the constant matrix (a n ^(h 3 ). For any A _L e 3 , the solution of (2.4) 
does not depend on the variable y 3 and then 

(W* - A • y) = and div (E n (fc 3 )VW*) = in 3t'(R 2 ). 

This equation and (2.4) imply that, for any A, [i _L e 3 , 

{a n )^h)X ■ ii = (£ n (/i)W n A ) • n = (£ n (/i 3 )W n A } • Jl = (a n )^h 3 )X ■ Ji. 

Hence, by Remark 2.2, (a n ) (h 3 ) converges to the 2x2 matrix <7* involved in (2.9). A two- 
dimensional perturbation formula in [ 7] gives the influence of the magnetic field h 3 on cr* : 

?. :=a*(ha) = a 2 + o^ 1 /^!) + h 3 PiJ, (2.11) 

where a® is a locally Lipschitz function defined on (0, oo) 2 , and for any a,\,a.i > 0, cr®(ai,ct2) is the 
transversal homogenized conductivity in the absence of a magnetic field. The independence of the 
microstructure of the variable x 3 allows us to obtain an explicit expression of a*{h) in terms of the 
transversal homogenized conductivity <r£ in the absence of a magnetic field. 

Remark 2.4. In the case where the high conducting phase is a set of circular fibres, it was proved 
in [17] that 0^(0:1,0:2) = 0:1/2 an d the limit a*(h) in (2.9) reduces to 

a*(h) = aih + a 2 + fi\o~^fk- Jh e 3 <g> e 3 + PiS(h). 

Now, let us proceed with the proof of Proposition 2.1. 

Proof of Proposition 2.1. Thanks to Remarks 2.2 and 2.3, there exists a 3 x 3 matrix a*{h) such 
that, up to a subsequence, we have the convergence of constant matrices 

(<r n ),(h) ><r*(h):=r* P *), (2.12) 
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where a* is given by (2.11) and where the constants q*,p* G R 2 , a* G R have to be deter- 
mined. To this end, we divide the proof into two steps. We first apply Grabovsky and Milton's 
method [21, 22, 23] to link (a n )^ to a more simple problem. Then, we study the asymptotic behav- 
ior of the different coefficients of this new problem. 

First step: A stable transformation under homogenization. For a fixed n G N*, following Grabovsky 
and Milton [23, 21], we consider two vectors Po,mQo,n G R 2 and the transformation 

~ / 5Vi P'n \ 

a' n := Yi n o n {h) n n = , (2.13) 



ot 



where 



and 



h 0\ ~ ^ / h Po,r. 

A J' U i 



n„ := ( T J , n» := L 1 > ( 214 ) 



p' = J a iPo,n- PiJh ^ mQ\Q n , q> = ) a i q °> n + Jh intt\Q n , {2 15) 

I &2,nP0,n 

Let us choose the parameters £>o,n and <7o,n in such a way that p' n and are constant. To that 
aim, po,n anci Qo,n have to satisfy the identities 

J 5lP0,n - /Si J/l = a 2 ,nP0,n ~ fa, n Jh, 

1 o f fqo,n + Pi Jh = aj n q 0yn + /3 2 , n Jh, 

which implies that 

P0,n = (/32,n - Pl)(&2,n ~ ^i)" 1 Jh and g ,n = {fo,n ~ Pl)(^l ~ &2,n)~ T Jh. 

The new matrix- valued function a' n defined by (2.13) is periodic and can be rewritten 

VzGf), a' n (x)=K(-) where K'-=\Z (2-16) 

V £ ™/ \q n a nJ 

Moreover, by (2.13), the coefficient a' n in (2.16) has the following explicit expression: 

y l,n = «1 + &lP0,n ■ Q0,n + Pi (po.n ~ ?0,w) " Jh, 
•^"'l.JrW+^n 1 ^ where ^ _ ( 2 - 17 ) 

a 2,n = Q 2,n + &2,nP0,n ' Q0,n + fo,n (P0,n ~ qO,n) ' Jh- 

Let us now study the homogenization of a' n . Define {(T n ) as in the formula (2.3). The conduc- 
tivity Y,' n does not depend on the variable y%. On the one hand, as in Remark 2.3, if A _L e$, the 
solution of the problem (2.5), with the conductivity T/ n , does not depend on the variable ys 
and VW^ = (VW^,0) T . Hence, since q' n is a constant, and by Remark 2.3, 

K)*A = ((S„W„ A ),(9; • W n A )) T = ({a n )~\,q' n ■ ~\f '. (2.18) 

On the other hand, it is clear that, for A = e^, W% 3 (y) = 1/3 satisfies (2.5) with the conductivity ~E' n . 
Hence, since p' n is a constant, we have 

K)*e 3 = ((p' n ,a' n ) T ) = {Pn, (Of- (2-19) 
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Then, by (2.17), (2.18), (2.19) and since \uj t , 

K)* : 



9 n , the matrix (o - ^)^ has the form 



\<Tn 



/T 

n 



Pn 
K), 



where 



"1 + <?lP0,n ■ 00,n + Pi (P0,n ~ qo,n) ' Jh 

«2,n + &2,nP0,n ■ 00,n + P2,n (P0,n ~ 0.0, n) ' Jh + o(l). 



(2.20) 



(2.21) 



Second step: Application of the theory of exact relations and asymptotic behavior of {o^j . By (1.3) 
and since the volume fraction 6 n converges to 0, we have 



Po. 



n — w n\ 



i{P2,n-Pl)(0n&2,n-a\)) ^ Jh > P 2 0~ 2 1 Jh 

00,n = n (/3 2 ,n ~ Pl){0n(o-1 ~ 0- 2 ,n)j Jh > -fh°2 Jh - 

Then, by (1.3), (2.15), (2.14), (2.21) and (2.22) we obtain the convergences 

' Pn -7*. P'* ■= [-Plh + 02*1*?] Jh, 



(2.22) 



n— too 

On > 0'*-=\ Plh ~ $2o{o 2 

n— ¥oo L 



Jh, 



(2.23) 



n, 



«) 



-»• n 



^ \ai - Pi a 2 X GiG 2 x Jh ■ Jh + 2p 2 Pi a 2 l Jh ■ Jh 

I 2 p 2 o 2 x Jh 



i=l 



and II r 



-»• n 



o 



(2.24) 



K p 2 h?Ja 2 l 1, 

Since the matrix transformation (2.13) preserves the H-limit in the periodic case (see, for in- 
stance, [21, 22, 29]), we have 

= n n (<r n ) t (h) n n . (2.25) 

Passing to the limit, as n goes to infinity, in relation (2.25), using (2.12), (2.20), (2.23)-(2.24), we 
obtain 

( a ; T ? ;)=na,n. (2.26) 

\0 * a J 

Inverting the identity (2.26) and taking into account (2.23) and (2.24) leads to (2.10). The proof of 
Proposition 2.1 is completed. □ 

Now let us turn to the non periodic case. 

3 The non periodic case 

In this section, we study the homogenization of the problem (1.4) without any periodicity as- 
sumption. The conductivity a n (h) is defined by 



a n (h) := a n I 3 + p n £(h) where \ n « , n « 

[ Pn := in\fi n "l + l n n P2,n- 

Consider the covering of M 3 by the squares Q\ defined by 

VfceZ 3 , Q k n = e n (Y + k). 



(3.1) 



(3.2) 



We assume that the conductivity coefficient a n defined by (3.1) satisfies, for any k € Z 3 , n € N*, the 
following conditions: 
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(z) the weighted Poincare-Wirtinger inequality 
VveH l (Q k n ), [ a. 



dx < Cn I a n \Vv\ z dx, 

Ik 



where c n is a sequence of positive constants satisfying 

c n > 0; 

(ii) there exists a positive constant c such that, for any fc£Z 3 and n £ N*, 

a n < c. 



(3.3) 



(3.4) 



fa, 



(3.5) 



Remark 3.1. Note that, in the periodic case, the hypothesis (2.6)-(2.7) is a rescaling of (3.3)-(3.4) 
which, similarly to the periodic case, prevents from the appearance of non local effects in the limit 
problem. Assumption (3.5) ensures that the microstructure does not concentrate on a lower dimen- 
sion subset through the homogenization process since it implies that (see in the proof of Lemma 3.1) 



r 1 



d€L°°(n) weakly-* in M(Q). 
In the periodic case, (3.5) is clearly satisfied since 

a n dx = ||a n || L i(y) < c, 



(3.6) 



where a n is defined by (2.2) and 8 = 1. 
We have the following result: 

Theorem 3.1. Assume that (1.3), (3.2)-(3.5) are satisfied. Then, there exist a matrix-valued func- 
tion cj*(/i) and a subsequence of n, still denoted by n, such that the solution u n of the problem (1.4) 
converges weakly in Hq[Q) to the solution u of the conductivity problem 



• div (o"*(/i)Vu) = / in Q, 
u = on dO,. 



Moreover, the effective conductivity cr*(/i) in (3.7) is given by 



(3.7) 



cr*{h) :- 



aT a* 



(3.8) 



where cr* is the H(A4(£l) 2 )- limit of o~ n (h) in the sense of Definition 1.1, 8 G L oc {Q) is the weak-* 
limit of ^~ 1 lo n and 



P* 



[Pih + / 02^ 1 (5*-?i)] T Jh, 
k a* = oli + 8ct2 + /3| 1 (<?i + 0&2 — a^a^Jh ■ Jh, 



(3.9) 



and, for any i = 1,2, 



Ui :- 



oti -/3i h 3 
J3ih 3 at 
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Remark 3.2. The shape (3.2) of is purely technical and can be generalized into any subset the 
diameter of which is of order e n . 

Remark 3.3. Since fl n has a columnar structure, ln n does not depend on the variable X3. Therefore, 

°n lt n n 9 e weakly-* in M(n). (3.10) 

Hence, as in Remark 2.3, it was proved in [17] that there exists a function cr£ defined on (0, oo) 2 and 
a subsequence of n, such that, for any ol\, a<i > 0, ^i,02 G 

H(M(n) 2 ) _ 

a n (h) = a n (h 3 ) — >■ a*(h) = cr°(a!i, a 2 + a 2 /3 2 /i 3 ) + h^fiiJ. 

We obtain, once again, an explicit expression of o*{h) in terms of the homogenized perturbed con- 
ductivity in the (x\, X2)-plane, in the absence of a magnetic field. 

A crucial ingredient of the proof of Theorem 3.1 is the following compactness result: 

Lemma 3.1. Let a n be the sequence defined by (3.1) such that (1.3) and (3.3)-(3.5) hold true. 
Consider two sequences £ n G L 1 ($7) and v n G H l (£l) satisfying 

£ n — £ weakly-* in JV[(VL) and v n — v weakly in H 1 ^). (3-11) 

We assume that 

[ a-^nl 2 dx+ [ a n \Vv n \ 2 dx < c. (3.12) 
Then, £ G L 2 (Q) and we have the convergence, in the sense of distributions 

inv n — ^ & m (3.13) 

Remark 3.4. Note that Lemma 3.1 is false when the conditions (3.3) and (3.4) do not hold. This 
can be seen by considering the classical model example of non local effects in conduction due to 
Fenchenko-Khruslov [ ] and presented, for instance, in [13, 10]. For the reader's convenience, we 
give the main steps of the counterexample. Let Q := (—5, ^) 3 and fi n be the i-periodic lattice of 
thin vertical cylinders of radius n~ 1 e~ n2 . Let a n be the conductivity defined by (3.1) with a\ := 1 
and «2,n := t 1 "" 1 e 2n which satisfies (1.3) and (3.5). For a fixed / in L 2 (f2), let u n be the solution, 
in Hq(U), of the equation 

-div(a n Vu n ) = / in^(fi). 
For R G (0, |), let V n be the y-periodic function defined on M 3 by 



V„(y) := 



2 

if r < e~ n (region of high conductivity) , 

1 if r>£ 



An easy computation shows that the sequences £ n := a„V« n • e 3 and v n {x) := V n (nx) satisfy the 
assumption (3.12) and that v n weakly converges to the constant function 1 in H l (Q). Moreover, 
Briane and Tchou [ ] proved that 

. du n du dv ii-i 

in = a n - * f := h ^— weakly-* m M(il), 

0x3 0x3 0x3 
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where the weak limit u of u n in H^(Q) and the weak-* limit v of 
measures satisfy the coupled system 



u n in the sense of Radon 



-Au + 27T (u — v) 
d 2 v 



+ v — u 



u 



v(x', 0) = v(x', 1) 



/ in 0, 

in U, 

on 90, 

if x'g(-±,±) 2 . 



(3.14) 



Then, if / is non zero, u and v are two different functions. Therefore, the convergence (3.13) does 
not hold true since, by the strong convergence, up to a subsequence, of v n to 1 in L 2 (f2) and the 
weak convergence of 1q\o Vit n to Vn in L 2 (0), we have 



1 



dx* 



^ / /- ^n dv . „,._. 

— ^xl = - + — m^'C7. 
0x3 0x3 0x3 



Substituting the expression of v, in terms of u, in the first equation of (3.14) leads to a non local term 
in the equation satisfied by u. The Poincare-Wirtinger control (3.3)-(3.4) is fundamental to avoid 
such effects. In this example, (3.4) is false since (see [12] for more details) the optimal constant 
in (3.3) satisfies c n > c > 0. 

Proof of Lemma 3.1. On the one hand, by (1.3), the sequence a n is bounded in L l (VL) and 
then, up to a subsequence, weakly-* converges to some a G Moreover, the Radon measure a 

belongs to L°°(f2). Indeed, let <p G ^b(^) an d denote again by tp its extension to R 3 by setting tp = 
on R 3 \ f2. There exists a finite subset /„ of Z 3 such that 

n c |J Qt 

where Q„ is defined by (3.2). As <p is a uniformly continuous function, we have 

/ a n (p dx = V" / a n ip dx = S2 ip(e n k) / a n dx + o(l). (3.15) 



kel, 



By (3.5), we have 



^2\ip(e n k)\ ! a n dx<cJ2 \Qn\Wnk)\ = c\\<p\\ L i {Q) + o(l) . (3.16) 

fce/n Q ™ keln 



The weak-* convergence of a n to a, combined with (3.15) and (3.16) yields 



tp a(dx^ 



< cIMIl 1 ^), 



which implies that the measure a is absolutely continuous with respect to the Lebesgue measure 
and a G L°°(il). From (1.3) and the convergence of a n to a, we have 

On^Qn = {GnOL2,n)~ 1 {a n ~ ailn\nj — 1 9 ■= aj 1 (a - ai) G L°°(0) weakly-* in .M(ft), 
and then a = a± + #0:2. 

On the other hand, by the Cauchy-Schwarz inequality combined with (3.11), (3.5) and the con- 
vergence of a n to a\ + o<2 9, we have, for any ip G ^(Cl), 



(p£{dx) 



lim 

n— >oo 



£n</> dx 



< lim sup / a n dx j a n tp z dx 

n—>oo JQ 



< c / (ai + 9a2)p 2 dx < c\\a± + 9a2\\ 
Jn 



(3.17) 
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which implies that the limit measure £ of £ n in (3.11) is actually an element of L 2 (Q). 

We now prove the convergence (3.13). Let tp E and let I n be a finite subset of Z 3 such that 

supp ip C |J C fi, 

fceJn 

where supp 99 is the support of 99. For any w E (S2), define It; 6 " the piecewise constant function 
associated with the partition {Qn) k£l as follows: 



w 



In order to study the convergence, in the sense of distributions, of (£ n v n — £v) to 0, we rewrite it as 

inVn-iv = i n {v n -V^ n )+i n (v^ n -V £n )+i n V £n - £v (3.18) 

N' V ^ 

:=Pn :=?n :=r ™ 

and estimate each term of the identity (3.18) separately. 

Convergence of the term p n in (3.18). Thanks to the Cauchy-Schwarz inequality, we have 

2 

in{v n -v^ n )Lp dx < ( 



V 



"n^nl 2 dx, / 

<iMi a -EL^i a ^E/ c 

kein Qn keln Q 

< c n||^||L/ "n 1 !^! 2 dx / a n |Vf n | 2 dx 



dx 



(3.19) 



dx 



where the last inequality is a consequence of (3.3). Finally, the inequality (3.19) combined with (3.12) 
and the convergence (3.4) yield 



in{v n -U^ n )(f dx 



< c^/c^ ► 0. 



(3.20) 



Convergence of the term q n in (3.18). By the Cauchy-Schwarz inequality and (3.5), we have 



£n(vn* n -v 6 ")^ dx 



< 



E 



Q'n 



£,n <P dx 



(v n — v) dx 



< \\<P\\ 



V] \Qn\~\ I «n dxj J a n l \&, 



Q'n 



2 dx / |f ra — v\ dx 



< 



E 



a n dx, an 1 !^! 2 dxJ / (v„-i;) 2 dx 
fc \/ ./ofc 1/ ./ofc 



,rf ./qs iri Jot 



dx 



kel n Qn kel n 
-lie |2 j„ /" /„. „.\2 



<c a n |£ n | dx / (u„ - v) dx 
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which yields, by (3.12), 



Jo. 



tp dx 



< c \\v n — v\ 



L 2 (n)- 



(3.21) 



Since v n converges weakly to v in Hq(TI), by Rellich's theorem, up to a subsequence, v n converges 
strongly to v in L 2 (n). Hence, (3.21) implies that 



tn{v^ £n -v 6 ")^ dx ► 0. 



(3.22) 



Convergence of the term r n in (3.18). Consider, for any S > 0, an approximation ip$ G "^(O) of v 
for the L 2 (fJ) norm, i.e., 

\\v-il>s\\v>(a)=o(S). (3.23) 

The term r n in (3.18) writes 



(3.24) 



On the one hand, since tpg £n converges uniformly, as n goes to infinity, to ipg E ^.(O), the conver- 
gence (3.11) of £ n to £ implies that the second term and the third term in the right hand side of the 
equality (3.24) converge to in Moreover, by the Cauchy-Schwarz inequality and the fact 

that £ € L 2 (f2), we have 



— v)(p dx 



< IMIooII^Il 2 ^)!!^ - v\\ L 2( n) . 



On the other hand, following (3.21), we have 



< C || C^lloo ||^ 



L 2 (Q)- 



Hence, by (3.23)-(3.25), we have 



lim sup 



[inVn" ~ (v)<p dx 



< c\\ips - v\\l 2 (q) = o{8), 



(3.25) 



(3.26) 



for arbitrary 5 > 0. 

Finally, putting together (3.18), (3.20), (3.22) and (3.26), we obtain that 



lim sup 



{inVn ~ &)(p dx 



0(5), 



which concludes the proof of Lemma 3.1. □ 
In the sequel we apply Lemma 3.1 to sequences £ n of vector-valued functions in L 1 (fi) 2 or L 1 (il) 3 . 

Proof of Theorem 3.1. Thanks to the equi-coerciveness a n > a\I^, the solution u n of the prob- 
lem (1.4) satisfies the convergence, up to a subsequence, 

u n — u weakly in Hq(TL), (3.27) 

for some u in iJg(Q). Moreover, putting u n as a test function in the equation (1.4), we obtain that 

/ a n |Vu„| 2 dx= a n Vu n ■ Vu n dx = (f,Un) H -i(n),Hi(n) < c - ( 3 - 28 ) 
J r2 J r2 
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Since ot 2 l l &2,n = I3 lajift.n^C')) by (1.3) the sequence |a 2 n cr 2,ri| is bounded. Then, as the se- 
quence a n is bounded in L 1 (0), the Cauchy-Schwarz inequality and (3.28) give 



2 

2 



a n Vu n \ dx < c / a n dx / a n \Vu n \ dx < c. 



Hence, we have the convergence of the current a n Vu n , up to a subsequence, 

a n Vu n — ^ £ weakly-* in M(Q) 3 , (3.29) 
for some £0 G .M(O) 3 . Moreover, by the boundedness of |«2~rj0"2,n| an< ^ (3.28), we have 

/ a~ l \a n Vu n \ 2 dx < c a n \Vu n \ 2 dx < c. 
Jn Jn 

Then, by Lemma 3.1 applied to £ n := <7 n Vit n , the measure £0 is actually an element of L 2 (S7) 3 . 

The rest of the proof, which is divided into three steps, is devoted to the determination of the 
form of the limit current £o- To that end, we use a method in the spirit of H-convergence of Murat- 
Tartar which is based on the cylindrical nature of the microstructure and the compactness result of 
Lemma 3.1 for sequences only bounded in L 2 (U; a n l ^ 2 dx). In the first two steps, we compute the 
components £0 • e\ and £0 ■ e 2 by combining Lemma 3.1 with a corrector function associated with 
the transversal conductivity a n , the existence of which is ensured by the two-dimensional results 
of [14, 17]. Since the corrector function considered in the previous steps is independent of the 
variable X3, the component £0 • 63 needs a different approach. This is the object of the last step. 

First step: Building a corrector. Thanks to Remark 3.3, up to a subsequence, a n H(Ai(Cl) 2 )- 
converges to some coercive matrix-valued function cr*. Then, the sequence i/(.A/f (f2) 2 )-converges 
to aj (see Theorem 2.1 of [14]). Let A G R 3 with A J_ e 3 . For A = (Ai, A 2 ) T G R 2 , let be the 
solution of 

div (5 T V^) = div(ajA) in 

v n = A • x on <9f2. 

By Definition 1.1, we have the convergences 

vfl — 1 A • x = A • x weakly in i? 1 (f2), 
ctJVdJ — ^ <x T A weakly-* in M(£l) 2 . 

Setting, for if!!, v^(x\, x 2 , X3) = i>^(xi,x 2 ), we have the convergences 



(3.30) 



(3.31) 



v 



A 



A-x weakly in H (fi), 



(3.32) 



[ «t t V^ — ^ aJX weakly-* in M(tt) 2 , 
and the energy inequality, as in (3.28), 

/ a n |V^| 2 dx = / Sy,V«* • Vv\ dx < a (3.33) 
Jn Jn 

Let tp G &>{£l). By (3.29) and since converges weakly to A • x in H 1 ^), putting v^cp as a test 
function in (1.4) yields 

/ cr n Vn„ • V (y*tp) dx ► (/, tp A • x) H -i (n) H i (n) = / £ • V(y?A • x) dx. (3.34) 

Jn n ^°° Jn 
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Since a n and do not depend on the variable X3, we have the identity 

a n Vu n ■ Vv^ = a n Vu n ■ Vv^ - d 3 (P n Vv* ■ Jh u n ) . (3.35) 
Then, by (3.35), an integration by parts gives 

/ a n Vu n -V(v^ip) dx = / a n Vu n ■ Vip dx (3.36) 
Jn Jn 



+ [ p n Vv* -Jhu n ^- dx (3.37) 
Jn ox 3 

+ [ a^Vv* ■ Vu n if dx. (3.38) 
Jn 



Step 2: Estimates of the terms in (3.36)-(3.38). The convergence of these terms are consequences of 
Lemma 3.1 and the generalized two-dimensional div-curl lemma in a high-contrast context of [13]. 

Convergence of the term on the right hand side of (3.36). On the one hand, by the boundedness 
of a n ^o~ n and (3.28), we have the inequality 

/ a-VnVMn-V^| 2 dx<c||V^||^ f a n \Vu n \ 2 <c. 
Jn Jn 

On the other hand, the convergence (3.29), the inequality (3.33), and the convergence (3.32) of u„ 
to A • x, show that the sequences £ n := a n \7u n and v n : = u„ satisfy the assumptions of Lemma 3.1. 
Hence, we obtain 

/ a n Vu n ■ V<p dx > / £0 • V</5 A • x dx. (3.39) 

Jn n ^°° Jn 

Convergence of the term in (3.37). We first compute the limit of f3 n \7v x in the sense of Radon 
measures. We have the identity 



— Ji -r2\f2„^i v ' 

where 

<Ti := I 

,Pih 3 ai J \P2, n h3 a 2 , n 

By (1.3), we have 



(3.40) 



o-i == „ , and a 2 , n := I ', . (3.41) 



~_T ~_T ~ I a 2 — $2^3 1 

/V^n > ft^a , where <j 2 := . (3.42) 

Combining this convergence with the ones in (3.31), we obtain that 

PrN^n — 1 0lA + fljcy 1 [5? - °i] A weakly-* in .M(ft) 2 . (3.43) 
By the boundedness of a" 1 ^, (1.3) and (3.33), we have 

/ a-^nV^I 2 dx<c [ a„|V^| 2 < c. 

J!) 

This inequality together with (3.28), (3.43) and the weak convergence (3.27) of u n to u in Hq(Q) 
show that the sequences ^ n '. — /3 n ^v^ and v n \ — u n satisfy the assumptions of Lemma 3.1. Then, 

f /3 n Vu* -Jhu n ^- dx ► f {Pil 2 + (3 2 [a* - a x ] a 2 l ) Jh ■ A u dx. (3.44) 

Jn 0x3 J Q dx 3 
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Convergence of the term in (3.38). Integrating by parts in (3.38), we obtain that 



aJVu^ • Vu n tp dx 



a^Vv^ ■ Vp u n dx + / o~n • V (u n p) dx. 



Tf/A 



(3.45) 



On the one hand, the boundedness of a n l a n and (3.33) yields 



a„ l \aEVv^\ 2 < c I a n \VvCV < c 



"n I n v 



,A|2 



By the second convergence of (3.32), the weak convergence (3.27) of u n to u in Hq(Q) and (3.28), 
the sequences £ n := a^Vv^ and w ra := u n satisfy the assumptions of Lemma 3.1. Hence, 



o~ n Vv n ■ Vp u n dx > I a* A • Vp u dx. 



n 



(3.46) 



in 

On the other hand, since a^Vv^ does not depend on the variable X3, the second term on the right 
hand side of (3.45) can be rewritten under the form 

•1 



5£Vi£ • V (u n p>) dx = /_ a n r V<(x') • V 



n 



(u n (p) (x',x 3 ) dx 3 







dx' 



(3.47) 



where x' = (xi, X2). 

In order to study the asymptotic behavior of (3.47), we apply a two-dimensional div-curl lemma 
of [17] which is an extension to the non symmetric case of [13]. Set, for any x' € CI, 



r) n {x') := a^Vv^(x') and v n (x') : = / {u n ip) (x',x 3 ) dx 3 . 

Jo 

Due to the convergences (3.31) and (3.27), we have 

r] n — 1 <tJA weakly-* in M.(£l) 2 



(3.48) 



V n X' 



f (x') := / (uc^) (x', x 3 ) dx 3 weakly in i^f 1 (^2). 



(3.49) 



The convergences (1.3) and (3.10), the definition (3.30) of the corrector (3.31)-(3.33) and (3.49) 
imply that the sequences rj n and v n defined in (3.48) satisfy the assumptions of Lemma 2.1 in [17]. 
Then, 



3£v^(x')-V 



(u n p) (x',x 3 ) dx 3 



aj\- V 



(wp) (x',x 3 ) dx 3 



in (3.50) 



Let -0 S f^(f2) such that i/> = 1 on the projection of the support of 99 on the (xi, X2) -plane. Taking ^ 
as a test function in (3.50), we obtain 



'_a£v^(x')-v 



(n n v?) (x',x 3 ) dx 3 



dx' 



> / <rjA • V (up) . 



Finally, this convergence combined with (3.45), (3.46) and (3.47) gives 

?JVw^ • Vm„ dx > I aj\ ■ Vu p dx. 

: n ->°° Jn 

Putting together (3.39), (3.44) and (3.51) with the equality (3.36)-(3.38), we obtain that 
0"nVn n • V(v^p) dx > / £0 • V95 A • x dx + / <t*Vu ■ A <p dx 

+ [ {Pih + h [?* - 3=i] S^ 1 ) fk-\u dx. 
7n <9x 3 



(3.51) 
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Since <t* depends only on the variable {x\,X2), this convergence, an integration by parts and (3.34) 
give 



£o • A <p dx 



, . Qui - ■ 

5*Vu - - — (Pxh + fa [5* - 3ql a? 1 ) Jh 
0x3 



A <p dx. 



(3.52) 



Finally, since the equation (3.52) holds for any ip £ ^(O) and any A _L e^, we obtain the first two 
components of £0 

£ = a*Vu - — {fal 2 + fa - ai] o^ 1 ) Jh. (3.53) 
Step 3: Computation of £0 " e 3- By (3.29), we have the convergence 



Oi n q V /3nVu n • Jh = a n Vu n • e 3 — >• £0 • e 3 weakly-* in M(tt). 

0x3 



(3.54) 



We first study the asymptotic behaviour of a n d3U n (which also gives the limit of fa d3U n due 
to the fact that, by virtue of (1.3), «2,n an d fa )U are of the same order). On the one hand, 
since 6 n = |fi| _1 |r2 n |, by the convergence (1.3), we have 



Q— !|0 



7ra«2,ra VnO!2,n 



< C. 



On the other hand, by (3.28), the weak convergence (3.27) of u n to u in Hq(Q) and (3.6), the 
sequences £ ra := 8~ 1 tQ n and v n := u n satisfy, once again, the assumptions of Lemma 3.1. Hence, 

On^QnUn — ^ OU in 

Moreover, since ln n does not depend on the variable x 3 , we have 



du r 



0x3 



(3.55) 



Finally, thanks to (3.55) and (1.3), we obtain the convergences, in the sense of Radon measures, 
( 

n>- = rvi II r-,V r-, 

dx 3 



9u n du n - du n 

Oin ~ = ttlln\n n ~ 1" \PnOt2,nW n l(l„ 



8x3 



9X3 

du n 



(ai + Va 2 )-z — , 

8x3 



0x3 



(3.56) 



8x3 ^-"^ 8x3 ^-"-;»>'« -»» dx3 
Now, in order to obtain the limit of the term fa\/u n in (3.54), which similarly to (3.40)-(3.41), writes 



/3nVli n = fat n \ nn Vu n + /32,n<7 2 ,n ^n^U n - Iq^Oi Vu„ 



(3.57) 



it remains to estimate a n Vu n . Since £o is the limit of the current a n \7u n (3.29) and since 

. . ^ du ~ • — - 

VA _L e 3 , a n Vu n • A = a n Vu n ■ X- fa -^-Jh ■ A, 

cte 3 

the equality (3.53) gives 

a n Vu n — fa —^-Jh — ^ <7*Vii — - — (fah + fa [<?* — ffi] cr^ 1 ) weakly-* in .M(O) 2 . 
0x3 0x3 



Then, combining this convergence with (3.56), we have 

a n Vu n — ^ a*Vu + fa — — ([ai + 9a 2 — or*] a? 1 ) Jh weakly-* in M(£lY 

0x3 



(3.58) 
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Finally, passing to the limit in (3.57), taking into account (3.58) and (3.42), we obtain the conver- 
gence, in the sense of Radon measures, 

, , , , 

f3 n Vu n — >■ [Pih + fav^fi* ~ Vn + ^2 + ^ 2 ~~ v*)°2 J h. 

Putting together (3.54), (3.56) and (3.59) yields 

£o-e 3 = [fhh + h^2 l {^* ~ ^i)] T J h - ^ u 



(3.59) 



+ 



du 
dx 3 ' 



(3.60) 



(at + a 2 9) + /3|cr 2 1 (a 1 + 9a 2 - a*)<r 2 1 Jh ■ Jh 

Finally, since the current a n Vu n weakly-* converges to £o in (3.29), we have the limit equation 

-div(£ ) = /, 

where, by (3.53) and (3.60), 

Co = (Co, Co • e 3 ) T = cr*(/t)V« 
which yields to the expression (3.8)-(3.9) of a*(h). Theorem 3.1 is proved. □ 

4 Two examples 

In this section we present two examples where the perturbation formulas for the effective con- 
ductivities of non periodic high-contrast columnar composites are fully explicitly computed. 

4.1 Circular fibres with variable radius 

Let p be a continuous function on Q depending only on the variable x' = (x%, x 2 ) satisfying 

3ci,C2>0, ci < p(x') < c 2 , Vx = (x 1 , 23) G Q, (4.1) 

and let r n be a sequence of positive numbers converging to 0, as n goes to infinity. We assume, 
without loss of generality, that 



/„ 

Jn 



p dx = 1. 



(4.2) 



Define, for any k £ Z 3 , the sequence (r n ^)neN* by 



r n ,k '■= r n \/p(e n k). 

We consider the case where Q n is the set of circular fibres = {y G Y \ y\ + y\ < r\ k } (see 

Figure 4.1) 



f2 n = f2n |J e n (uj nyk + k). 
Note that the fibres io n ^ do not have the same radius. 



(4.3) 



2r, 



n,k 




Figure 4.1: The cross section of the non periodic microstructure 
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We have the following result: 



Proposition 4.1. Let Ct n be the sequence of subsets of CI defined by (4.3) and a n (h) be the associated 
conductivity in the problem (1.4). Assume that 



-> 0. 



(4.4) 



Then, there exist a matrix-valued function a*(h) and a subsequence of n, still denoted by n, such that 
the solution u n of the problem (1.4) converges weakly in Hq(CI) to the solution u of the conductivity 
problem 

- div (o\t(/i)Vu) = / in CI, 
u = on dCl, 



where <7*(/t) is given, for any x = (x',x 3 ) E ; fry 



<7*(/t)(a;) = ai/ 3 + p(x ) 



/a^ + a 2 /3 2 2 |/i| 2 



«2 + /^l 



e 3 ®e 3 + /?i<r(/i). 



(4.5) 



Remark 4.1. We can easily check that the homogenized conductivity cr*(h) of the two dimensional 
microstructure of the Figure 4.1 is given by 

a*(h) = a\h + P\h 3 J = a\{h). 

This leads to the simple form (4.5) of a*(h). 

Proof of Proposition 4.1 In order to apply Theorem 3.1, we need to check that the condi- 
tions (3.3)-(3.5) hold true. On the one hand, the Poincare-Wirtinger inequality combined with (1.3) 
imply the existence, for any £ Z 3 , of a sequence of positive constants c n ^ such that 









2 




/ a n 


v — r v 


JQ 



a n \Vv\ dx, 



(4.6) 



where Q n = £ n (Y + k). Using estimates derived in [12], one can show that the best constant in the 
weighted Poincare-Wirtinger inequality (4.6) satisfies 



V/c G 1? , Vn G N*, < < ce 2 In (r„ y/p(e n k)) 
for some positive constant c. Therefore, by (4.1) and (4.4), we have 

< c n k < ce n \ hxr n \ + o(l) > 0, 

n— >oo 

uniformly with respect to k G Z 3 . Conditions (3.3) and (3.4) of Theorem 3.1 are satisfied. On the 
other hand, by the definition of Cl n and (4.2), we have the following estimate for the volume fraction 



I Cl n | 

\Cl\ n— yoo 



7757 J2 4 7rr nP( £ nfc) ~ irr 2 n f pdx = vrr 
\il\ n-toc / 



which, by (4.1), implies that for any n G N* and k G Z 3 , 

-f a n dx = ai(l -7rr^p(e n k)) + a 2 , n ^r n p(e n k) < c + c9 n a 2 , n < c. 

jQ k n 

Then, condition (3.5) of Theorem 3.1 is satisfied. Theorem 3.1 and Remark 4.1 ensure the existence 
of an effective conductivity a*(h) which, after an easy computation, writes 



a*(h)(x) = a\h + 0(x') 



, ( a\ + a 2 fil\h\ 2 



e 3 ®e 3 + pi£(h) Vx = (x',x 3 ) G Cl x (0,1), (4.7) 
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where 8 is the weak limit of 8 n 1 lfi n . The function 9 in (4.7) coincides with p. Indeed, since p is 
continuous, we obtain, for any tp G ^o(JT) extended to M 3 by setting ip = on M 3 \ 17, 



/ °n dx = —2^2[ V dx + = —2 ^^nP^M^) + o(l) 

which implies that 



/ ^n 1 lnn^ d ^= / P^dx + 0(1). 

in in 



Finally O' 1 !^ converges weakly-* to p in and, then, 8 = p. This concludes the proof of 

Proposition 4.1. □ 

4.2 Thin squared grids 

In this section, we consider the case of a columnar composite the cross section of which is a highly 
conducting grid surrounded by another conducting medium (see Figure 4.2). Let t n be a positive 
sequence converging to as n goes to infinity. Let p be a continuous function on f2, depending only 
on the variable x' = (x\,X2) and satisfying 



3ci,C2 > 0, Cl < p(x') < C2, Vx = (x',xs) G 17. 

We assume, without loss of generality, that 

p dx = 1. 



Define, for any k in Z 3 , the sequence (i n 



/ 



fcJneN" 



by 



Let Q, n be the set of non periodically distributed squared fibres 

Q n = Q, n (J e n (uj n:k + fc) where := {y G V | max(|yi|, |y 2 |) > 5 - Vfc} • 

Note that the case p = 1 leads to a periodic distribution of the squared fibres in f2. 



(4i 



(4.9) 



(4.10) 





(a) Periodic case 



(b) Non periodic case 



Figure 4.2: The cross section of the structure 



We have the following result: 
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Proposition 4.2. LetQ n be the sequence of subsets ofVt defined by (4.10) anda n (h) be the associated 
conductivity in the problem (1.4). Assume that 



->• «2 > and 4t n /?2' 



(4.11) 



Then, there exist a matrix-valued function o~*(h) and a subsequence ofn, still denoted by n, such that 
the solution u n of the problem (1.4) converges weakly in Hq(£1) to the solution u of the conductivity 
problem 

- div (<j* (h) Vw) = / in O, 
u = on d£l, 



where a*{h) is given by 



and, for any (x',x^) £ ; 



(T*(h) 



a* (/i) 



<7*(/i) = ai + p(x ) 



p* 



2a 2 



jh + p ( x ')^h, 



2a 2 



/3i + /j(x' 



J/i + p(x ) 



2a 2 



h, 



(4.12) 



(4.13) 



(4.14) 



a* = ai + /j(x') a 2 + p(x') r-^~(^i + 
In formula (4.14), p = 1 corresponds to the periodic case. 

Proof of Proposition 4.2. Let us first consider the periodic case. In order to apply Proposi- 
tion 2.1, we need to check that (2.6) and (2.7) are satisfied. To this end, consider V £ ^(Y) such 
that (V) = 0. Define, for any n E N*, the subsets K l n , i = 1,2, 3, 4, of Y by 

K\ := [—g, 2] x [2 ~~ 2] X [ — 2> 2] ' ^« ' = [ _ 2' 2] X i~2 + — 2] X [ — 2' l] ' 
:= [2 ~~ 2] X [ _ 2' 2] X [ _ 2> 2] ' "^n := [ _ 2 + ~~ 2] X [ _ 2' 2] X [ _ 2' 2] ' 

For instance, the projection of K\, in the (yi, y 2 )-plane, is the shaded zone in Figure 4.3. 




?72 

4- 




Figure 4.3: The period cell of the cross section of the microstructure 
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By the definition (2.2) of a n , we have 

( a n V 2 dy < [ a n V 2 dy + V f a n V 2 dy 

I— 1 * * 

< ax f 4 V 2 dy + a 2 , n y] [ V 2 dy 

1=1 

< a\ [ V 2 dy + c^V f V 2 dy. 

Since (V) = 0, this inequality and the Poincare-Wirtinger inequality in i7 1 (y), yield 

f a n V 2 dy< ai [ \VV\ 2 dy + o 2 ,„ V f V 2 dy. 
Jy Jy i=1 Jxi 



(4.15) 



We now estimate the second term of the right hand side of this inequality. On the one hand, since K\ 
is convex, the Poincare-Wirtinger constant in H l {K^) is bounded from above by the diameter of K\ 
divided by ir [ ] and, therefore 



V dy < 2 
Kl JKl 



V--t V dy 

Kl 



+ 2\Kh 



V dy 



< c 



I / |W| 2 dy+\Kl\ i V dy 
\JkI Jki 



Kl 
2 N 



(4.16) 



On the other hand, for any — ^<r, s,t<^,we have 

V(s) - V{r) = [ V'(t) dt, where V(t) := 



1 i 

2 / 2 



1 11 

2 2 



V(yi,t,y 3 ) dyidy 3 . 

Integrating the first equality in (4.17) with respect to s £ [5 — t n , i] and r 6 [—5, i], we have 



(4.17) 



V dy — -f V dy 

Kl Jy 



< 



2 / 2 / 2 



2 ln 2 2 



|y'(t)| dt < 





dV 




dy 2 



dy, 



which, since (V) = 0, implies 


that 






i V dy 






JKl 





(4.19) 



< j |W| dy < \\W\\ L 2 (Y)3 . (4.18) 
Then, combining (4.16) and (4.18) with the boundedness (4.11) of 0:2,71 = in«2,n, we obtain that 

«2,n V 2 dy < c ^a 2 , n | W| 2 dy + ||W \\ 2 L 2 {Y) A < c j^a n \VV\ 2 dy. 
Similarly to (4.19), we have, for i = 2,3,4, 

«2,n f V 2 dy < c [ a n \VV\ 2 dy. 
Jki Jy 

Finally, (4.15), (4.19) and (4.20) imply 

a n V 2 dy<cj a n |W| 2 dy. (4.21) 



(4.20) 



Y 
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By a density argument, (4.21) is satisfied for any V £ H l (Y) with (V) = 0. Since e n converges to 0, 
the hypotheses (2.6) and (2.7) of Proposition 2.1 are satisfied. Then, there exists a homogenized 
matrix which is given in terms of the transversal effective conductivity a* of the microstructure 
of Figure 4.2a. It remains to determine 5"*. Since one can choose the cross-like shape of the Fig- 
ure 4.3 as the period cell of the transversal microstructure of the heterogeneous medium occupying f2, 
Proposition 3.2 of [17] ensures that 



Ct\ + 



2a 2 



I 2 + Pih 3 J, 



(4.22) 



and formula (4.13)-(4.14) of cr*(h) is a consequence of (2.9)-(2.10) where <7* is given by (4.22). The 
periodic case is then proved. 

The existence of cr*(h) in the non periodic case is a consequence of Theorem 3.1. Indeed, for 
any k £ 1? and n 6 N*, a rescaling of (4.21) gives 



Jo] 







2 




V — -T V 


dx < ce 2 n \ 


'Q k n 




JQ 



\Vv\ dx, 



and, by (4.9), 



■f a n dx = ai(l -4t„ jfc (l - f H)fc )) + 4a 2 , n t n ,k( 1 ~ f n,fc) < c + ct n a 2 , n < c. 

The assumptions of Theorem 3.1 are satisfied. Then, there exist a matrix- valued function a*{h) and 
a subsequence of n, still denoted by n, such that the solution u n of the problem (1.4) converges 
weakly in Hq(Q) to the solution u of the conductivity problem (4.12). In view of the formulas (3.8) 
and (3.9), the expression of cr*(h) becomes explicit as soon as a* and 9 are identified. 

On the one hand, it is easy to check, similarly to the proof of Proposition 4.1, that p is the 
weak-* limit, in the sense of Radon measures, of the sequence 6~ l \^i n . Then, the function 6 in 
Theorem 3.1 turns out to be p. On the other hand, by Remark 3.3, in order to compute cr*, one 
has to determine a®(ai, a 2 ), which is the H(Ai(Q) 2 )-lim.it, in the sense of Definition 1.1, of the 
conductivity cr ra (0), in the absence of a magnetic field, given by, for any x' € £1, 



5 n (0) := 



a\h in \ ^n, 
ct2,nh in Sin- 



Due to the local nature [13] of the .ff (.M(f2) 2 )-convergence, it is sufficient to compute cr®(cti,a 2 ) lo- 
cally in £1. To that aim, consider x' G £1 and e > small enough such that the closed disk D{x' , e) C O. 
Since p is continuous and by (4.8), we have 

sup p{z). (4.23) 



< c h£ (x') := mf p{z) < p(x') < c 2 , e (x') := 

z£D{x',e) 









D(,v' 


LJ 


'/////A 


V//////// 


' t 






H 


LJ 












(a) Bound from below (b) Bound from above 

Figure 4.4: Bounds from below and above of <J n (0) 
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For % = 1,2, let rt n be the subset of D(x',e) defined by (see Figure 4.4) 

Qif = D(x',e)n (J e n (k + {y e (-3, \) 2 | max(|yi|, \y 2 \) > I - <v(a/) t n }) , (4.24) 

fcgZ 2 

and let <t^ be the periodic conductivity defined on D(x',e) by 

Sj. := J"i*» in gf^)\^ (425) 

By the definitions (4.23) and (4.25), we have for any z £ D(x',e), the inequalities 

^ e (z)<a n mz)<a^{z). (4.26) 

For the rest of the proof, we need the following result which is a consequence of the two-dimensional 
div-curl lemma, in a high contrast context, of [13]: 

Lemma 4.1. Let D be a bounded domain ofM 2 and, fori = 1,2, consider an equi-coercive sequence 
of symmetric matrix-valued functions A l n £ L°°(D) 2x2 bounded in L 1 (D) 2x2 which H(M.(D) 2 )- con- 
verges to A\ in the sense of Definition 1.1. We assume that 

VnGN*, A n <A 2 n a.e.mD. (4.27) 

Then, we have the inequality 

A] < A\ a.e. m D. 
Proof of Lemma 4.1 Let A G K 2 . Consider, for i = 1, 2, the solution Vn l of 

J div(^V^' i ) = div(4A) in D, 
{ Vn' 1 = A • x on dD. 

By Definition 1.1, we have the convergences, for i = 1,2, 

{Vn % — ^ A-x weakly in H^(D), 
Affivn'* — >■ A* A weakly-* in M(D) 2 . 

On the one hand, by (4.27), we have the inequality, almost everywhere in D 

2 A n Vv x n l ■ V^< 2 - A n Vv^ ■ Vv x n l < A n Vv x n > 2 ■ Vv x n 2 < A 2 n Vv x n 2 ■ Vv^' 2 . (4.28) 

On the other hand, applying, for i,j = 1, 2, the two-dimensional div-curl lemma of [13] (Theorem 2.1) 
to £ ra := A n Vvn l and v n := Vn 3 we have the convergences, in the sense of distributions, 

Vi, j = 1, 2, £ n ■ Vv n = A n Vv^ ■ Vv^ A\X ■ A in &(D). (4.29) 

Finally, combining (4.28) and (4.29), we obtain 

2Al\ - X — AlX ■ A < A 2 X ■ A in &'(D), 

which concludes the proof of Lemma 4.1. □ 

Since, for i = 1,2, a % n is an equi-coercive sequence of periodic matrix-valued functions bounded 
in L 1 (D(x', £■)), a % n H(Ai(D) 2 )-convevges to a constant matrix al' £ . Then, applying Lemma 4.1 
with D = D(x',e) and (4.26), we have 

< o°(a 1 ,a 2 )(z) < a 2 / a.e. z e ^(x',e). (4.30) 
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Moreover, due to the definition (4.24) of Qlf and the convergence (4.11), we have 
\ttlf \ \ fl\' 1 a 2 , n = 4:t n a 2 ,nCi, e (x') + o(l) >- a 2 Ci t e(x') > 0. 

n— ¥oo 

Then, substituting a 2 Ci tE (x') for a 2 in (4.22) in the absence of a magnetic field (i.e., = 0), we 
obtain, for i = 1,2, 

a 1 + c i Jx')^)l 2 . (4.31) 



2 

By (4.23) and (4.31), taking the limit, as e goes to 0, in the inequalities (4.30), we obtain, for any 
Lebesgue point x' of cr® (011,012) in f2, 

<T°(a 1 ,a 2 ) = (cti + p{x') ^) h- 
Therefore, by Remark 3.3, we have 

5* = U + p(x') d+MH) h + fiihzJ. (4.32) 

Finally, we apply the formula (3.8)-(3.9) for <r*(/t) in Theorem 3.1, with a* given by (4.32), to 
obtain (4.13)-(4.14). This concludes the proof of Proposition 4.2. □ 

Acknowledgements: The authors wish to thank M. Briane for his valuable comments and suggestions. 
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